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The proton matrix element of the isovector-scalar density, (p\uu — dd\p)/2M p , 
• is calculated by evaluating the nucleon current correlation function in an 

external isovector-scalar field using the QCD sum-rule method. In addition 
Qh| to the usual chiral and gluon condensates of the QCD vacuum, the response 

Q_|| of the chiral condensates to the external isovector field enters the calculation. 

The latter is determined by two independent methods. One relates it to 
the difference between the up and down quark chiral condensates and the 
other uses the chiral perturbation theory. To first order in the quark mass 
difference 5m = m^ — m u , the non-electromagnetic part of the neutron-proton 
mass difference is given by the product 5m(p\uu — dd\p) /2M p ; the resulting 
value is in reasonable agreement with the experimental value. 
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I. INTRODUCTION 



Understanding the properties of the nucleon is of obvious importance in hadron physics. 
A variety of nucleon matrix elements of bilinear quark operators have been evaluated in the 
past using various approaches. The QCD sum-rule method originally proposed by Shifman, 
Vainstein, and Zakharov [Ij and its extension by Ioffe and Smilga 0] for external field prob- 
lems, provide a tractable framework for the study of these nucleon matrix elements. These 
include the matrix element of electromagnetic current to determine the magnetic moments 
0-^], the matrix element of the axial vector current to find the renormalization of nucleon 
axial coupling constant []5]|| , the matrix element of the quark part of the energy momentum 
tensor, which gives the momentum fraction carried by the up and down quarks in deep 
inelastic scattering [0,0, and the matrix element of isoscalar-scalar current for evaluating 
the nucleon Sigma term ||. In the present work, we calculate the matrix element of the 
isovector-scalar density, (p\uu — dd\p) /2M p , following the same external field approach. 

The appearance of the external field leads to specific new features in QCD sum rules 
which distinguish them from those in the absence of the external field. Thus the phenomeno- 
logical representation of the correlation function written in terms of physical intermediate 
states contains a double pole at the nucleon mass whose residue contains the matrix element 
of interest. In addition there are single pole terms which arise from the transition matrix 
element between the ground state nucleon and excited states. These later contributions are 
not exponentially damped after Borel transformation relative to the double pole term and 
should be retained in a consistent analysis of the sum rules. In the theoretical side of the 
sum rules expressed in terms of an operator product expansion (OPE) the external field 
contributes in two different ways-by directly coupling to the quark fields in the nucleon 
current and by polarizing the QCD vacuum. 

Since for our problem the external field is a Lorentz scalar, non-scalar correlators cannot 
be induced in the QCD vacuum. However the external field does modify the quark and gluon 
condensates already present in the QCD vacuum. It turns out that for the problem under 
study the most important one is the response of the up and down quark chiral condensates 
to the external field which can be described by a susceptibility x- This can be determined 
by writing a spectral representation, which can be evaluated for example by using chiral 
perturbation theory. Alternatively, since in the real world up and down quark masses are 
different, 5m = ma — m u itself can be regarded as an external isovector scalar field and % 
can be related to the isospin breaking in quark condensates 7 = ((0\dd — uu\0)) / (0\uu\0) 
and is given to first order in 5m by x = —7 /5m. At the current level of accuracy where 
these two determinations can be done, they are mutually consistent. 

In Sec. [Tl] we derive the sum rules for the nucleon current correlation function in an 
external isovector-scalar field and describe the analysis of these sum rules, which leads us to 
determine (p\uu — dd\p) /2M p . We also study its dependence on the x value. 

In Sec. |TJ we show that an evaluation of the matrix element (p\uu — dd\p) /2M P enables 
us to compute the non-electromagnetic part of the neutron-proton mass difference. Since 
the 1970 's it has been recognized that the empirical mass difference arises from two sources. 



One is purely electromagnetic and yields a contribution of —0.76 ± 0.30 MeV [10[] to the 
neutron-proton mass difference. The second source is the difference between up and down 
quark masses, which is of the order of the quark masses themselves. This difference leads to 
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larger contribution which overrides the electromagnetic contribution making neutron heavier 
than proton. To first order in 5m this later contribution is given by the product of 5m and 
(p\uu-dd\p)/2M p . 

Recently several authors [ll]-|14| have applied the QCD sum rules for the nucleon mass 
to extract the neutron-proton mass difference by including the up and down quark masses 
in the mass sum rules for proton and neutron and the isospin breaking in the condensates. 
The relation of these works to ours is described in Sec. [TV], where we also briefly comment 
on the Nolen-Schiffer anomaly. 



II. SUM-RULE CALCULATION 

The procedure that we use for calculating the matrix element (p\uu — dd\p) /2M P follows 
the same pattern as used in Refs. @-|9[]. We first couple the quarks to an external isovector- 
scalar field Sy(x) which is described by adding a term AC to the usual QCD Lagrangian 

AC = -S v (x)[u{x)u(x) - d{x)d{x)} . (2.1) 

We can take Sv{x) to be a constant Sy since we are only interested in the zero momentum 
transfer matrix element. The correlation functions of the nucleon current in QCD vacuum 
in the presence of Sy will be computed. The term linear in Sy gives the matrix element of 
interest . 



A. QCD sum rules for (p\uu — dd\p) 



Consider the correlation function of the proton interpolating field in the presence of a 
constant external isovector-scalar field S v which can be taken to be arbitrarily small 

(2.2) 



(2.3) 



U(S v ,q) eez / d*xe^(0\T[r) p (x)rj p (0)]) Sv 



where r/ p is the proton interpolating field introduced in Ref. flT5 

v%(x) Cjf.u^x) 757 M 4(z) , 



x 



^abc 



where u a (x) and d c (x) stand for the up and down quark fields, a, b and c are the color indices, 
and C = —C T is the charge conjugation matrix. Since Sy is a scalar, Lorentz covariance 
and parity allow one to decompose n(SV, q) into two distinct structures 



U(S v ,q)^U\S v ,q 2 )+W(S v ,q 2 )<i . 



(2.4) 



To the first order in the external field S V} the two invariant functions can be written as 



W(S v ,q 2 ) = Ul(q 2 ) + S v U\(q 2 



(2.5) 



for {i = l,q}, where I1q are the invariant functions in the absence of the external field which 
give rise to the mass sum rules discussed extensively in Refs. [p!5|-[T8| . We are concerned here 
with the linear response to the external field given by U\(q 2 ). 



3 



To derive a QCD sum rule, one first carries out an OPE, which will express U\(q 2 ) 
in terms of various vacuum correlators, and then matches it to an expansion in terms of 
physical intermediate states. Now the external field contributes to U(S v ,q) in two ways: it 
couples directly to the quark fields in the propagating nucleon current and it also polarizes 
the QCD vacuum. The chiral condensates of up and down quarks change as follows 

(uu) Sv = (uu) - xSv(uu) , (2.6) 
(dd) Sv = (dd} + X S Y (dd) 0j (2.7) 

where (O) = (0|O|0). Using Eq. fl2.1|) it is easy to see that x is related to the correlation 
function 

X(uu) = l -J d 4 x(T{w(x)w(x) - d{x)d{x),u{0)u{0) - 3(0)d(0)})„ . (2.8) 

Similarly the mixed quark-gluon condensates change as follows 

(g s lla ■ Qu) Sv = (g s ucr ■ Qu) - XmSv(g s ucr ■ Qu) , (2.9) 

{g s da ■ Qd) Sv = (g s da ■ Qd) + XmSv(g s d<7 ■ Qd) , (2.10) 

where a ■ Q = a^ v Q^ u with Q^ v the gluon field tensor and x m , the susceptibility corre- 
sponding to the quark-gluon mixed condensate, can also be expressed in terms of a spectral 
representation. 

Since isospin is a good symmetry for hadron matrix elements we have assumed that the 
response of the up and down quarks is the same, apart from sign, and that we can disregard 
the vchange in the gluon condensate ((a s /n)G 2 ) due to the external isovector field. 

To calculate the Wilson coefficients of the OPE, we need the coordinate-space quark 
propagators in the presence of the external field and the vacuum condensates. To first order 
in the external field S v , the propagators in the fixed-point gauge |T9|-f21~f take the form 

(T[<(^(0)]) 5v = JL**1 - - ^(uu)^ 

+^ ab X S v {uu) Q 5 i3 + S ab -^(uu) [#] 
_ ^192^ 0r ' ~ Sah XmS v —(9sU(T ■ Qu) Q 5ij 

— ^ (G»MT b ~ 2 IM U + ^tk + ■■■, (2.H) 



32vr 2 v " X 2 
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-^ ab xS Y (dd) 5 l3 - 5 ab ^(dd) ]ft 

2 2 

^192 {9s2a ' Gd)A3 + ^^^Ssda ■ gd)A, 
-^f^{9s~da-Gd) [t) i3 

~— 2 {G,MT b - 2 {^ y + tTfto + ■■■■ (2.12) 



Here we are interested in terms linear in S v . We can then disregard the current quark masses 
because they make negligible contributions. 

The correlation function n(S v , q) can be computed using the Eqs. (|2.11| ) and (|2.12|) 
above. We have computed the contributions corresponding to the diagrams listed in Fig. |l|. 
The results of our calculations for the invariant functions 11^ and II} are 

nf(g 2 ) = -^(dd) \n(-q 2 ) + - (2(g s ua ■ Qu) a - (g s da ■ Qd)^j — 



+ ^(uu) (g s ua ■ Qu)J- + ^-(uu) (g s ua ■ Qu) \ , (2.13) 
o g 4 o g 4 



nj(g 2 ) = ^i^fH^) + ^- 2 (dd) q 2 H-q 2 ) - 2 - (3(uu) (dd} - 2{uu) 2 ) i • (2.14) 

We now turn to the phenomenological side of the sum rules, which is obtained by ex- 
panding n(SV, q) in terms of physical hadronic intermediate states. There are three types 
of contributions. Firstly the matrix element of interest is contained in the term 

(0\v P \p)(p\uu-dd\p)(p\v P \0) , (2-15) 

where the current rj p creates a proton that interacts with the external field S v and is then 
annihilated by the current r} p . Defining 

(0\V P \P) = \v , (2.16) 

where X p denotes the coupling between rj „(0) |0) and the physical proton state and v is the 
usual Dirac spinor (vv = 2M p ), and introducing the notation 

_ {p\uu - dd\p) 
H= 2M P ' (2 ' 17) 

one can write the contribution of Eq. fl2.15|) to Tl\(q 2 ) as 
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It is seen that the above term has a double pole at the nucleoli mass. 

The external field S v can also cause transition between the proton and an excited state 
which can have either positive or negative parity relative to the proton. When the relative 
parity is positive, the contribution can be written as 



XpXp* 



+ M V 



H* 



i-M* 



Q 



Ml " q 2 - M* 



(2.19) 



where now H* is the transition matrix element between \p) and \p*). This term has a simple 
pole at the proton mass as well as at the mass M* of the excited state. It is easy to see that 
after a Borel transformation this contribution is not exponentially damped as compared to 
the double pole contribution Eq. Q2.18p . Therefore, one has to retain these simple pole terms 
in the analysis of the sum rules through the introduction of a phenomenological parameter 
to be determined along with the diagonal matrix element H (see Refs. |2|-||). The third 
type of contributions comes from transitions involving only the excited states. These are of 
course exponentially damped after a Borel transformation and can be approximated in the 
usual manner [|i5|-|ITj|j , by equating them to the perturbative contributions starting from an 
effective threshold. 

Equating the OPE results Eqs. ( J2.14|) and ( J2.13| ) and the physical intermediate state 
expansion discussed above and applying the Borel transformation [|TJ, we obtain the following 
sum rules 

^- 2 (qq)oE L-^ + \ X (m)lL A " - ^(M^ ■ Gq) L~ 8/9 - ^(qq)o(g s qv ■ Gq) L~ 2/27 



Xm /_ 



6M 2 



(qq)o(gsqa-Gq) L- 2 / 27 



2 M P 

P JP H + q 



-M?/M 2 



[2.20) 



M 6 



EnL 



-8/9 



+ 



X 

4vr 2 



(qq)oM*E 1 --(qq) 2 



r, M 2 



-Mp/M 2 



(2.21) 



where Ax and A q are the phenomenological parameters that represent the sum over 
the contributions from all off-diagonal transitions between the proton and the excited 



states. Here we have defined E = 1 - e~ s °/ M , E x = 1 - e- s °' M [j^ + lj and 

£"2 = 1 — e _s °/ A/2 + -f^s + Ijj which account for the sum of the contributions involv- 
ing excited states only, where Sq is an effective continuum threshold. In Eqs. ( |2.20|) and 
( |2.21| ), we have ignored the isospin breaking in the vacuum condensates. We have also 
taken into account the anomalous dimension of the various operators through the factor 
L = \n(M 2 / Aq CD ) / ln(;U 2 /AQ CD ) |])|l5j. We take the renormalization scale /1 and the QCD 
scale parameter Aqcd to be 500 MeV and 150 MeV. 
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B. Estimate of x 



It is clear from Eq. ( J2.8[ ) that x is determined once the isovector-scalar two point function 
is known. The latter has been studied by Gasser and Leutwyler using chiral perturbation 
theory to one loop |2"2fl . They found 

i [ d 4 x(0\T{u(x)u(x) - d(x)d(x),u(0)u(0) - d(0)d(0)}|0) = 8 ( — ^ — \ h 3 . (2.22) 
J ■ \m u + m d ) 



Since two pions cannot form an isovector-scalar, the | 7T7t) intermediate state does not con- 
tribute in Eq. ( |2.22|) . The other possible pseudoscalar two particle states are \KK) and \r]ii) 



which means that an extension to SU(3) flavor symmetry is necessary. This has been done 



by Gasser and Leutwyler in Ref. p3| . Using Eq. (11.6) of Ref. p3[] , and (ss) /(qq) = 0.8, 
we get h 3 ~ —0.003. Combining Eqs. ( |2.8p and ( |2.22| ), we obtain 

2 

777 

(m u + m d )f* 

where we have used (m u + m d )(qq) = —m%f%, and the experimental values m n = 138 MeV, 
f w = 93 MeV, and a median value (qq) = — (240MeV) 3 which corresponds to m u + m d = 
11.8 MeV. Alternatively, x can be determined as follows. The terms proportional to the 
current quark masses in the QCD Lagrangian can be written as 

Anass — — 77i (uu + dd) + -5m(uu — dd) — m s ss — ■ ■ ■ , (2.24) 

where rh = \{m u + m d ) and the ellipses denote the terms due to heavier quarks. Treating 
5m(uu — dd) as a source term one obtains using Eq. ( |2.8|) 

x(uu) = -—{uu-dd) . (2.25) 
dom 

On the other hand, one can expand (uu — dd) and - d j^(uu — dd) into the Taylor Series in 
5m. Using (uu — dd) \s m =o = 0, we find 

(uu - dd) = x$m{uu) + 0[(5m) 2 ] , (2.26) 

which implies 

x 5m = - 7 + 0[(5m) 2 ] . (2.27) 

Therefore, to the lowest order in 5m, the susceptibility x is determined by the ratio of the 
isospin breaking parameters 7 and 5m. The value of 7 has been estimated previously in 
various approaches [p3} -p9|, with results ranging from — 1 x 10~ 2 to —3 x 10~ 3 . Gasser and 



Leutwyler ||2^] have determined the ratio 5m/(m u + m d ) = 0.28 ±0.03. Since we have used a 
median value of 11.8 MeV for the sum of the up and down quark masses we adopt a median 
value for 5m = 3.3 MeV. Here we consider the 7 values to be in the range —1 x 10~ 2 to 



-3 x 10 3 , which, upon using Eq. ( |2.27|) , corresponds to 



0.9 GeV" 1 < x < 3 GeV" 1 . (2.28) 

The susceptibility Xm can also be determined using a spectral representation for Eqs. ( |2.9| ) 
and ( |2.10|) . However, given the uncertainty in the value of x> we assume, in this work, 

Xm — A • 
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C. Sum-rule analysis 



Defining a = -4TC 2 (qq) , A 2 = 32vr 4 A 2 , m 2 , = (qg s a ■ Gq) /(qq) , A q = (2n) 4 A q , and 
Ai = (2tt) a Ai, we can rewrite the sum rules Eqs. ( |2.21| ) and fl2.20| ) as 



,M p 2 /M 2 



-M 4 aE L-^ 9 + 



AM 2 



M 



L 4/9 + £-8/9 _ ^ a 2 m 2 L -2/27 

6 3 



\\M V H + A q M 2 



(2.29) 



M 8 E 2 L- s / g - M 6 X aE 1 - ±M 2 a 2 



\;M 2 H + A\M' 



(2.30) 



To extract H = (p\uu — dd\p) /2M P , A q and A\ from the above sum rules, we use the 
experimental value for the proton mass M p and extract A 2 from the proton mass sum rules. 
Using a = 0.55 GeV 3 (m u + = 11.8 MeV) and ml = 0.8 GeV 2 , one finds from the best fit 
of the proton mass sum rules that X 2 = 2.1 GeV 6 corresponding to so = 2.3 GeV 2 @.We use 
only the first or the chiral-odd sum rule of Iofle which is more accurate than the second or 
the chiral-even sum rule. 

First consider the sum rule Eq. (|2.29|) , which is obtained from Ul(q 2 ). For definiteness 
we take for \ the value \ = 2.2 GeV -1 given in Eq. ( |2.23| ). In Fig. 0, we plot the individual 
terms in the LHS of Eq. (|2.29|) as well as their sum, as functions of M 2 in the interval 
0.8 < M 2 < 1.4 GeV 2 . It can be seen that the LHS follows a linear behavior in M 2 and we 
can match it with the RHS. To find the best values for the constant and the coefficient of the 
linear term in the RHS we follow the numerical optimization procedure used in Refs. [T7| , |3"0"| ] . 
We sample the sum rules in the fiducial region of M 2 , where the contributions from the 
highest- dimensional condensates included in the sum rule remain small and the continuum 
contribution is controllable. Here we choose 0.8 < M 2 < 1.4 GeV 2 as the optimization 
region, which is identified by Ioffe and Smilga @ as the fiducial region for the nucleon mass 
sum rules. To quantify the fit of the left- and right-hand sides, we use the logarithmic 
measure 



5{M 2 



In 



maximum{Li/S', RHS} 
minimum{LifS', RHS} 



(2.31) 



which is averaged over 100 points evenly spaced within the fiducial region of M 2 , where 
LHS and RHS denote the left- and right-hand sides of the sum rules respectively. The 
sum-rule predictions are obtained by minimizing 5. Using this procedure we obtain 



H ~ 0.54 



A ~ 
i\ q — 



0.29 GeV J 



(2.32) 



The RHS of Eq. ( [2.29|) with these optimized values is also plotted in Fig. 2. One can see that 
the LHS and RHS have a very good overlap. 

We now turn to the other sum rule Eq. ( |2.30| ), which is obtained from U\. In Fig. |3|, 
the individual terms in the LHS as well as their sum are shown as functions of M 2 for 
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X = 2.2 GeV -1 . The LHS is fairly linear in M 2 in the interval 0.8 < M 2 < 1.4GeV -2 . The 



RHS of Eq. ( 2.30 ), with the optimized values 

H ~ 0.01 , A x ~ -1.82 GeV 6 , (2.33) 

is also shown in Fig. || 

It is worth noting that in the sum rule Eq. (|2.29j) the double pole term is more important 
than the single pole term, which is also qualitatively evident from the near constancy of the 
LHS as a function of M 2 . In contrast, the single pole term in the sum rule Eq. (|2.30|) clearly 



dominates, leading us to suspect that it is not reliable to determine the double pole term in 
which we are interested. This is confirmed by the following. 

Let us consider relaxing our tacit assumption that the continuum threshold used in our 
external field sum rules should be the same as the one occurring in Ioffe's mass sum rules. 
If sq is varied from 2.3GeV 2 to 2GeV 2 or 2.6GeV 2 , the result for the matrix element H 
extracted from the sum rule Eq. (|2.30|) changes from 0.01 to —0.07 (s = 2GeV 2 ) or +0.08 
(sq = 2.6GeV 2 ) while the result from the sum rule Eq. ( |2.29| ) changes only from 0.54 to 0.52 
or 0.56. 

It is clear from the above analysis that the chiral-odd sum rule Eq. ( |2.29| ) is extremely 
stable, while the chiral-even sum rule Eq. ( |2.30| ) is not. So, in this paper we shall disregard 
the results based on the sum rule Eq. (|2.30|) and consider only the results from the stable 
sum rule Eq. (|2.29| ). 

The fact that one of the sum rules works well, while the other fails, is not peculiar to 
the problem under study. This pattern is seen also in a study of the isoscalar-scalar matrix 
element as well as in the sum rules for the matrix elements of electromagnet current and 
axial vector current [0-0]. As discussed extensively in Ref. 0, the different asymptotic 
behavior of various sum rules can be traced to the fact that even and odd parity states 
contribute with different sign and kinematical factors. If chiral symmetry is realized in the 
Wigner-Weyl mode at high energies, i.e., by parity doubling, it is possible to have either 
cancellation or reinforcement between excited state contributions. Irrespective of the exact 
manner in which these cancellations take place, it is clear that the sum rule in which the 
continuum contributions are weak is more reliable. This is the case for the chiral-odd sum 
rule Eq. ( |2.29| ), where the continuum factor Eq appears, to be contrasted with the chiral-even 
sum rule Eq. (|2.30|) , where the factor E 2 occurs. 



Let us now consider the effect of varying x > which as we have seen earlier is not precisely 
known. We find that the quality of the overlap of the two sides of the chiral-odd sum rule 
remains good (as measured by 5) as we change the value of \i an d gets better as x increases 
and worse as x decreases. We also find that the continuum contribution gets larger as x 
decreases and smaller as x increases. For x values in the range 1.4 GeV -1 < x < 3.0 GeV -1 , 
we find 

H = 0.32 - 0.76 . (2.34) 

For x < 1.4 GeV -1 , we find that the continuum contribution is larger than 50%. This 
implies that the sum rule is dominated by continuum and the predictions are not reliable 
for x values smaller than 1.4 GeV -1 . 
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III. THE NEUTRON-PROTON MASS DIFFERENCE 



In this section we consider the relation between the matrix element evaluated in the last 
section and the neutron-proton mass difference. First let us disregard electromagnetism and 
work to the first order in the quark mass difference m^ — m u . 

Consider the quark mass term in the QCD Hamiltonian density H. QCU , as given by 

'Hmass = rn u uu + m d dd + m s ss H 

= m(uu + dd) — -5m(uu — dd) + m s ss + • • • . (3-1) 

The isospin symmetry is explicitly broken by the term proportional to 5m. Using covariant 
normalization for the hadron state labeled by h and momentum k 

(k',h\k,h) = {2Txfk°5^ ] {k'- k) , (3.2) 



and regarding the 5m term in Eq. ( |3.1| ) as a small parameter, we can make a perturbation 
expansion and write the shift in the mass Mh of the hadron as [|K| 



5(M h ) = Z^m^^R . (3.3) 

The matrix element occurring in Eq. ( [3. 3D is to be computed at 5m = 0, which means isospin 
can be taken to be exact. We can then write for the difference between the neutron and 
proton mass to first order in 5m as 

(M„ - M r ), = sJj^m . p.4) 

The subscript q in the left-hand side denotes the fact that we are considering only the 
non-electromagnetic part of the mass difference. 

The effect of turning on electromagnetism is described by an effective Lagrangian 



d A yD{x-y)Tf{x) 3 ,{y) , (3.5) 



where D(x) = [i4n 2 (z 2 — ie)]^ 1 is the photon propagator and j^{x) is the electromagnetic 
current. To remove the divergence arising from electromagnetism one must add counter 
terms, and these, of course, depend on renormalization prescription. However, this depen- 
dence is extremely weak. For example the change in the up quark mass due to a change in 
the renormalization scale by a factor of two is less than 0.01 MeV. It is therefore meaningful 
to separate the contribution from the quark mass difference, from that due to electromag- 
netism. The latter is estimated to be |L0] 



(M n - M p ) elec = -0.76 ± 0.30 . (3.6) 
The experimental mass difference is 1.29 MeV, which then gives 
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(M„ - M p ) e q xp = 2.05 ± 0.30 . (3.7) 

In the last section we saw that uncertainty in our knowledge of x leads to a corresponding 
uncertainty in our determination of H. For the \ value obtained from chiral perturbation 
theory [see Eq. ( gjp ], we get (M n - M p ) q ~ 1.8 MeV. 



For x values in the range 2.15 GeV 1 < x < 2.80 GeV 1 , we find 

1.75MeV < (M n - M p ) g < 2.35MeV , (3.8) 

which is consistent with experimental data. Smaller and larger values of x outside the range 
considered above lead to correspondingly smaller and larger values for the neutron-proton 
mass difference. 



IV. DISCUSSION 

One of our main objectives in this paper has been to extract the proton matrix element 
H = (p\uu — dd\p) /2M p . We have seen that the chiral-odd sum rule Eq. ( |2.29| ) is reliable for 
determining this matrix element. The major limiting factor has been the uncertainty in the 
value of the susceptibility x- A more accurate evaluation of the two point function Eq. ( |2.8| ) 
should reduce the uncertainty in the value of x an d hence help to pin down the value of the 
matrix element H. 

We also saw that the non-electromagnetic part of the neutron-proton mass difference 
is essentially given by the matrix element H multiplied by the light quark mass difference 
5m. If we use a median value 5m = 3.3 MeV and H ~ 0.54 as obtained using a value 
of x — 2.2GeV _1 we get (M n — M p ) q ~ 1.8 MeV, which indeed has the right sign and 
magnitude. This suggest that our approach to extract H and the neutron-proton mass 
difference is reliable. However, since x an d 5m are not precisely known we cannot make a 
critical comparison with data at present. 

We now turn to a comparison of our method with those of earlier authors. The nucleon 



mass was originally extracted by Ioffe [H| with a combined use of both the chiral-odd and 
chiral-even mass sum rules. It is necessary to use these two sum rules together since one 
must eliminate the coupling constant A^. Belyaev and Ioffe [16[] extended this method to 



determine the mass splitting between hyperon and nucleon by treating the strange quark 
mass as a perturbation. In addition to the mass shift, My — M^, one must also take into 
account the change in the coupling constant A y — Xj^ and the change in the continuum 
threshold Sq. The authors of Refs. flTB|JT4| used the same procedure to determine the mass 



splittings within an isospin multiplet by treating rrid and m u as perturbation parameters. 
In Ref. fll[ the neutron-proton mass was extracted directly from the difference between 



the neutron and proton mass sum rules, but the continuum contributions were disregarded. 



In Ref. [13], Adami et.al. retained continuum corrections but regarded 7 as a parameter 



to be determined by a fit to all isospin splittings in the baryon octet. In Ref. |T2|, apart 



from the perturbation due to the quark masses, an attempt was made to incorporate the 
electromagnetic contribution also phenomenologically in the sum rules. 

The sum rules derived by us in Sec. |T| can also be derived directly from the mass sum 
rules. Writing m u = m — 5m/2, m d = m + 5m/2 and using x = —j/5m one can differentiate 



11 



Eqs. (16) and (17) of Ref. jl2] with respect to 8m. One can then identify our sum rules 



Eqs. ( ggg) and fl2~30l) with Eqs. (29) and (30) of Ref. 0. (An assumption about the mixed 
chiral condensate, equivalent to our assumption, Xm = was made in Ref. [0). This 
coincidence between the sum rules is not surprising, since the quark mass term in the QCD 
Lagrangian can also be regarded as a constant external scalar field. The double pole term 
in the RHS of our sum rules clearly arises from the simple pole term in the mass sum rules 
after the differentiation of the proton mass with respect to 5m. 

We note that the term proportional to the mixed quark-gluon condensate [the fourth 
term in the LHS of our Eq. ( ^.29|) ] has not been included in Ref. [JT3]]. We have seen in our 
analysis of the sum rules (see Fig. 2) that this term is numerically significant. There is also a 
minor discrepancy in the coefficient of the four-quark condensate term between Ref. |13|] and 



this work (or Ref. |12|]). This discrepancy comes from the fact that the authors of Ref. [13 



directly used the £ and H mass sum rules from Ref. ||16|| , where not all the quark mass terms 
were taken into account. 

It is now easy to discussed earlier in Sec. why Eq. (|2.29|) is a better sum 

rule as compared to the chiral-even sum rule Eq. ( |2.30| ). In the chiral-odd sum rule, the 
perturbation due to the finite quark mass does not affect the leading asymptotic behavior 
and hence when a differentiation with respect to 5m (alternately when the neutron and 
proton sum rule difference is taken) the leading asymptotic behavior reduces to a smaller 
power of the Borel Mass M 2 , which in turn means weaker dependence on the continuum 
contribution. On the other hand, in the chiral-even sum rule the introduction of quark mass 
leads to the term m^M 6 in the proton and m u M e in the neutron sum rule. Consequently in 
the difference the leading asymptotic behavior is now M 6 . In other words, the continuum 
contributions are enhanced. This feature is also clearly reflected in our analysis in Sec. ||. 
For the chiral-odd sum rule the double pole term residue was stable when s was varied 
while in the chiral-even sum rule the corresponding residue was unstable. Further, in the 
chiral-even sum rule the single pole term corresponding to transition between proton and 
the excited states dominated over the double pole term. We would also like to point out 



that inclusion of instanton contributions improves the chiral-even mass sum rule |31| . In 
the light of our observations above regarding the relation of our sum rules to the mass sum 
rules it suggests that such instanton contributions can also be significant in our sum rule 
Eq. ( ggpp . 



It is also worth emphasizing that in our work two stages are involved. We first calculate 
a hadronic matrix element H. This involves the susceptibility Xi which being essentially the 
ratio —'j /5m is relatively insensitive to errors in our knowledge of 5m. The quark mass part 
of the neutron-proton mass difference was obtained as the product H5m. 

As a final remark we shall comment on the Nolen-Schiffer anomaly |32|]. We saw in the 
last section that the empirical neutron-proton mass difference can be written as 

(M n - M p ) exp = -0.76 + 5mH . (4.1) 

Now, it is well known that the matrix element of the axial vector current is quenched [3"5fl 
inside the nuclear medium by about 30 percent. It may be reasonable to assume that the 
isovector-scalar matrix element (p\uu — dd\p) /2M p is also quenched in a similar fashion. 

Assuming then for example a 30 percent reduction in the value of H, it follows form 
Eq. ( (4.1| ) that the effective mass difference in the nuclear medium is (M n — M p )^f d ~ 
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0.49 MeV. Understanding the Nolen-Schiffer anomaly is then reduced to explain the quench- 
ing of (p\uu — dd\p) /2M p in nuclear medium. This can be handled either by traditional 
nuclear structure calculations or again by use of QCD sum rules as in the quenching of 



nucleon axial coupling [34,35 
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FIGURES 



FIG. 1. Diagrams for the calculation of the Wilson coefficients of the correlation function. The 
solid, wavy and dashed lines represent the quark, gluon, and external field, respectively. 



dependence of the left-hand side (solid curve) and right-hand 
of Eq. Q2.29 ), with x = 2.2 GeV™ 1 and the optimized values 
dd\p)/2M p = 0.54 and A q = 0.29 GeV 5 . The curves 1, 2, 3 and 4 correspond to the 



FIG. 2. Borel mass 
side (long-dashed curve) 



first, second, third and fourth terms in the LHS of Eq. (2.29). 



FIG. 3. Borel mass dependence of the left-hand side (solid curve) and right-hand 
side (long-dashed curve) of Eq. (2.30), with x = 2.2 GeV™ 1 and the optimized values 
(p\uu — dd\p) /2M P = 0.01 and Ai = —1.82 GeV 6 . The curves 1, 2 and 3 correspond to the 
first, second and third terms in the LHS of Eq. (2.30). 
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